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ABSTRACT 



The influence of inelastic hydrogen atom collisions on non-LTE spectral line formation has been, and remains to be, a significant 
source of uncertainty for stellar abundance analyses, due to the difficulty in obtaining accurate data for low-energy atomic colli- 
sions either experimentally or theoretically. For lack of a better alternative, the classical "Drawin formula" is often used. Over recent 
decades, our understanding of these collisions has improved markedly, predominantly through a number of detailed quantum me- 
chanical calculations. In this paper, the Drawin formula is compared with the quantum mechanical calculations both in terms of the 
underlying physics and the resulting rate coefficients. It is shown that the Drawin formula does not contain the essential physics behind 
direct excitation by H atom collisions, the important physical mechanism being quantum mechanical in character. Quantitatively, the 
Drawin formula compares poorly with the results of the available quantum mechanical calculations, usually significantly overestimat- 
ing the collision rates by amounts that vary markedly between transitions. 
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1. Introduction 

Solution of the non-local thermodynamic equilibrium (non-LTE) 
radiative transfer problem in stellar atmospheres requires de- 
tailed and complete knowledge of the radiative and collisional 
processes that affect the statistical equilibrium of a given atomic 
species of interest. A difficulty regarding the collisional pro- 
cesses is to determine which, among the almost endless possi- 
bilities in a stellar atmosphere, are important. In an ear ly study 
of the formation of the Na D lines in the solar spectrum. iPlaskettl 
d 1955b considered the two obvious candidates for the case of 
the solar atmosphere: inelastic collisions with electrons and hy- 
drogen atoms in their ground state. The importance of electron 
collisions in many environments is well known, arising from the 
fact that electrons are always the most abundant charged parti- 
cle and have a much higher thermal velocity than atoms, leading 
directly to a higher collision rate. 

Hydrogen atoms are the most abundant perturber by far; 
however, as pointed out by Plaskett, given that in thermal con- 
ditions atoms are much slower than e lectrons, simpl e consider- 
ations based on the Massey criterion dMassevll 19491) indicate H 
atom collisions are nearly adiabatic and thus should have much 
smaller inelastic cross sectio ns than electro n collisions at the rel- 
evant thermal velocities (see lLambertl 1 993[ for a detailed discus- 
sion). Plaskett considered experimental estimates for quenching 
of Na by He and methane, which suggested deexcitation cross 
sections ranging from non-detectable for He to a small value 
for methane. Assuming that quenching by H might be similar 
to that for methane, Plaskett estimated the inelastic H atom col- 
lision rate would be of the same order of magnitude as that for 
electrons, despite n(H)/n(e) ~ 10 4 . As indicated by Plaskett, if 
one made the more reasonable assumption that H would have 



a similar cross section to He, then the rate would probably be 
negligible. 

Based on this rough expectation, and probably also because 
of the almost total lack of any data for atom-atom collisions at 
low energy due to the increased difficulties in calculating and 
meas uring such cro ss sections compared to the case of electrons 
(e.g. iMassevi [l955h . most early non-LTE studies of cool stars 
during the 1960's and 1970's neglected H atom collisions. While 
some autho rs mentioned the possible need to include such pro- 
cesses (e.g. Ivan Regemorter|[l962tlAihav & Canfieldll 19691) . the 
first to make an attempt at i ncluding them quantitatively was the 
study of solar Na lines by Gehren] (| 19751) . That study used the 
classical formula from iDraw in (1969), though details on how 
these formula were extended from the case of collisions between 
like atoms to collisions between unlike atoms are not given. 
Gehren found the H atom collisions to be relatively unimportant, 
less than half of the electron collision rates. 

A later study of Li line formation by Steenboc k & Holw eger 
(1984) used a similar approach, though in this case the general- 
isation of the formula of Drawin to inelastic H atom collisions 
was now given. The resulting formula is often referred to in stel- 
lar astrophysics as the "Drawin formula". The key development 
of this study was that it considered the case of a low-metallicity 
halo dwarf star, where n(H)/n(e) ~ 10 6 . They found that if H 
atom collisions were described to order-of-magnitude accuracy 
by this formula, then indeed the large ratio n(H)/n(e) would be 
sufficient to offset the lower efficiency and frequency of H atom 
collisions relative to electrons, and inelastic collisions with H 
atoms would be the dominant collisional process in cool metal- 
poor dwarfs. Thus, the study of Steenbock and Holweger restim- 
ulated the idea that H atom collisions might be important in cool 
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stellar atmospheres, and since that study the question has been: 
does the Drawin formula really provide order-of-magnitude es- 
timates of these processes? 

Recent theoretical and experimental studies have put us in a 
position to answer this question. In this paper we will examine 
the question from two points of view. First, is there any reason- 
able expectation that the Drawin formula can provide reasonable 
results given our present understanding of the physics involved? 
This is discussed in Sect. [2] Second, how do the results from the 
Drawin formula compare with those from quantum mechanical 
calculations? This comparison is made in Sect. [3] Finally, the 
main conclusions are summarised. 



2. Physics of inelastic hydrogen collisions 

In this section we discuss the physics behind the classical 
Drawin formula. We then compare with quantum theory and the 
physics of inelastic H atom collisions as revealed by experiment 
and detailed quantum scattering calculations based on quantum- 
chemistry descriptions of the relevant quasi-molecules. 



2. 1 . The classical Drawin formula 

The Drawin formula is the result of a number of modifications 
and extensions of the clas sical formula for ionization of atoms 
by electron impact due to Thomson ( 1912). In Thomson's the- 
ory, the bound electron in the target atom is considered as a sta- 
tionary free classical electron. The Coulomb interaction of an 
incident classical electron with the atomic electron is then con- 
sidered and Thomson calculates the impact parameter for which 
the deflection of the incident electron results in an energy trans- 
fer between electrons corresponding precisely to the ionization 
potential. All collisions inside this impact parameter will have 
a larger energy transfer and thus lea d to ion i zation and so the 
cross section can be easily calculated. lElwertl (1 19521) introduced 
a factor to account for equivalent electrons in the outer shell of 
the atom, which is usually also applied. 

As pointed out by Rudge ( 1968, Sect.4), Thomson's theory 
is remarkable for being able to make a number of qualitative pre- 
dictions regarding ionization cross sections and their behaviour 
with impact energy. However, quantitativel y it provides o nly an 
order-of-magnitude estimate. For example, Seaton (1962) com- 
pared with experimental results for H, He, Ne and Ar, and found 
the classical theory overestimates by a factor of about 5 at low 
energy. Seaton explained this from the fact that in quantum me- 
chanical treatments the ionization probability rarely approaches 
unity even at small impact parameters, while the classical the- 
ory assumes the probability is unity for all impact parameters 
smaller than that leading to a classical energy transfer equal to 
the ionization potential. It was also observed that the comparison 
differed markedly for Ne, a difference that is explained by quan- 
tum mechanical properties of the atom, and cannot be explained 
by classical theory. 

In Fig. Q] we complement Seaton's comparisons with new 
comparisons with modem experimental data from the Queen's 
University of Belfast laboratory for the cases of electron impact 
ioniz ation of ground state hydrogen, oxygen , magnesium and 



iron dShah et all I 987H Thompson 'etd]ll995t iMcCallion etail 



1992: IShah etal 



19931 respectively). These experimental data 



are in good agreement with other experiments, cover a wide 
range of collision energies including those near threshold and 
the absolute uncertainty is expected to be less than 20% (see 
caption to Fig. [TJ. The general features of the comparison are 
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Fig. 1. Comparison of the classical Thomso cross sec- 

tion for ionization of H, O, M g and Fe atoms by electron impact 
with e x perimental result s of IShah et al.l (Il987j): [T hompso n et all 
dl995h : lMc"Calrion et alJd!992l) : IShah et alJdl993l) . respectively. 
The plot shows the ratio of the experimental values with respect 
to the theoretical cross section, plotted against the collision en- 
ergy E scaled by the ionization energy /. The errors in the ex- 
perimental values are in all cases typically 5% or less, with an 
additional uncertainty of around 7-14% in the absolute scale. 



very similar to those seen in Seaton's comparison, i.e. the classi- 
cal formula overestimates the cross sections by around a fac- 
tor of 5 at low energy. In this comparison, O is seen to be 
particularly discrepant, differing by a factor of 20 near thresh- 
old. Thus, quantitively, the Thomson formula is able to provide 
only a rough order-of-magnitude estimate of the ionization cross 
section. Given the difficulties involved in quantum calculations 
compared to the classical approach, considerable efforts were 
made to improve the classical models, predominantly through 
inclusion of the initial velocity of the atomic electron, the ef- 
fects of the target nucleus, and exchange effects between elec- 
trons. These efforts were met with mixed s uccess, and there is 
an extensive literat ure on the subject (see iBurg ess & Perciva 
1968; lRudgdll968l iBelv & van Regemorterlll970t lBatesll!97 
for some reviews). 



iDrawinl d!968l 119691) : iDrawin & Emardl (1 19731) extended 
the Thomson theory to the case of ionization and excitation 
in collisions between like atoms, which was later generalised 
by ISteenbock & Holwegerl \ 1984) to the case of H atom col- 
lisions with other atom s. The d erivat ion has been redone in 
considerable detail by iLambertl (Il993l) . correcting an error in 
the cross section expressions due to confusion of centre-of- 
mass a nd laboratory frames; an error t hat was pointed out ear- 
lier by Fleischmann & DehmeJ (Il972l) . It is worth noting that 
if this error was carried through one would expect consider- 
able errors in the rate coefficients due to the incorrect threshold. 
However, the equations for th e rate coefficients generalised by 
ISteenbock & Holwegerl d 19841) seem to be corr ect, with the fi- 
nal expression for the rate coefficient derived bv lLamberj (119931) 
differing from that of Steenbock and Holweger by an essentially 
negligible factor ma / (m\ + toh), of order unity. The resulting ex- 
pression for the rate coefficient for an inelastic H atom collision 
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with atom A, A; + H — > A u + H, where I and u denote lower and 
upper states respectively, as presented by Lambert is: 

/ \ ia z( E f] f m e (m A + m H ) (2kT\ U2 

{(rv) = I6na — //„- ■ — T(w„ ; ) (1) 

[EfJ (m H + tn e )m H \ n/j ) 

where //„ is the oscillator strength for the transition / — > u, is 
the ionization potential of H, = E^ - E^, w„i = E^/lcT, and 
Y is a function arising from the integration over the Maxwell 
velocity distribution, which for kT ~ E^, is: 

= exp(-w a/ )(l + — |. (2) 

\ W u lj 

All other quantities have their usual meaning. 

The main physical assumption in the generalisation of the 
Thomson theory to ionization in collisions between like atoms, 
is that for a given collision energy above threshold, the efficiency 
of the energy transfer in an atom-atom collision is the same as 
in the electron-atom collision. This amounts to neglect of the 
nucleus of the perturbing atom, and thus the model in essence 
is the classical interaction of the impact atom electron with the 
stationary free classical atomic electron, allowing for the in- 
creased mass of the perturber. Drawin introduces an additional 
mass fact or (mA/mn) without justification , wh ich has b e en dis - 
cussed by [Fleischmann & Dehmell d 19721) and iLambertl (11 993). 
While Drawin compared with experiments for ionization in col- 
lisions between light atoms and between molecules and found 
agreement within a factor of two, Fleischmann & Dehmel made 
comparisons and found only order-of-magnitude agreement be- 
coming worse for heavier atoms. In the case of collisions be- 
tween atoms of similar mass, Fleischmann & Dehmel found only 
agreement within a factor of 100. Comparisons wi th more recent 
exper imental data do not change this conclusion dKunc & Soonl 
[1991 . 

Classical descriptions of atomic processes are less well 
suited to excitation than to ionization, due to the gener- 
ally smaller energy transfer and the difficulty in c hoosing a 
classical final energy s tate band. In dDrawinl |l969t see also 
iDrawin & Emardlll973l) the ionization formula is extended to 
the case of excitation by analogy with the ionization formula; 
that is, the ionization threshold energy is replaced by the excita- 
tion threshold energy, and the oscillator strength is inserted into 
the formula. This amounts to a choice of the final energy band 
to be all energy transfers greater than the excitation threshold, 
modified by the oscillator strength. The inclusion of the oscil- 
lator strength is almost certainly motivated by analogy with the 
Bethe approximation for inelastic collisions with electrons (see 
below). This approximation is valid only for optically allowed 
transitions and at high collision energies where the momentum 
transfer is small. The oscillator strengt h is introduce d in some 
versions of the ionization formula (e.g. lDrawi r3 fl969l) including 
the one adopted by Steenbock and Ho lweger, but not in others 
(iDrawinll 1 968t IDrawin & Emardll 1 9731) . 

2.2. Quantum collision theory 

In certain cases, quantum scattering theory leads to a depen- 
dence of the excitation (and ionization) cross section on the 
(generalized) oscillator strength. This arises within the first 
Born approximation (e.g. [ Seatonll962tlLandau & Lifshitzll965(: 
Bran sden & Jo achain 20031), which is valid for high-energy col- 
lisions and is obtained within perturbation theory. For electron- 
atom collisions, the unperturbed wave functions are taken in the 



form of atomic initial/final wave functions <f>f {u) (r) (r = {r,,} de- 
noting a set of coordinates of atomic electrons labeled by a) 
multiplied by plane wave functions (2ny 3 ^ 2 exp(tk/(„)R) of the 
impact particle (R being the impact electron coordinate) and the 
perturbation operator is the interaction potential U between an 
impact electron and the target atom, which consists of Coulomb 
interactions with the nucleus C/ en and the atomic electrons U ee : 
U = U en + U ee . The direct scattering probability amplitude for 
excitation of the atom / — > u is 

A B] = -^<C'(r)l exp( iq R){/|^f(r)> , (3) 
Zn 

where q = k/ - k„ is the momentum transfer, and the integration 
is taken over both r and R. Integration over the impact elec- 
tron coordinate R, considering the Coulomb interactions and the 
orthogonality of the atomic wave functions, removes the inter- 
action with the nucleus and leaves only the interactions between 
the incident and atomic electrons 

A B1 = -4<C'(r)l J] exp(iqr a )|0* (r)) r . (4) 

The cross section is proportional to the square of the scattering 
amplitude in Eq.|4] and finally the c ross section can be shown to 
be proportional to the quantity (e.g. lBransden & Joach ain 2003) 

r hl = 2(£ "~ £/) K^'(r)l 2 exp(;qr ff )|^ t (r)) r | 2 , (5) 

q a 

which is called the generalized oscillator strength for reasons 
that will become apparent. For optically allowed transitions, ex- 
pansion of the exponential in Eq.[5]can be restricted to the dipole 
term leading to a relation between the excitation cross section 
and both the transition dipole moment and the corresponding 
optical oscillato r strength fi u ; this relation is known as the B ethe 
approximation (|Befhd Il930t ISeatonl Il962t lLandau & Lifshitz 
ll%5l: iBransden & Joachainll2003l) . and leads to 

ISeatonl d 19621) provides a discussion of the physical interpreta- 
tion of this relationship between the collisional and radiative pro- 
cess. For optically forbidden transitions, the higher-order expan- 
sion terms resulting from Eq. |5]must be considered thus yielding 
nonzero excitation/ionization cross sections. Thus, the first Born 
approximation forms a bridge between the quantum approach 
and the classical approach of Drawin, both being based on the 
Coulomb interactions of an impact charged particle with atomic 
electrons and resulting in relations between the excitation (ion- 
ization) cross section and the oscillator strength. 

The Born approximation is applicable for high-energy 
electron-atom and electron-molecule collisions, but is valid also 
for charged heavy particle (ion) - atom/molecule collisions. For 
low-energy collisions, which are the main interest of this paper, 
the Born approximation is no longer valid and gives rather poor 
results (further explaining the poor agreement between the clas- 
sical calculations and the experimental data at low energies, see 
Fig. 1). 

At low collision energies, quantum atomic collision the- 
ory uses a different approach, in particular, the standard adi- 
abatic (Born-Oppenheimer) approach, which is based on the 
separation of the nuclear kinetic energy operator and the elec- 
tronic fixed-nuclei Hamiltonian in the total Hamiltonian, fi- 
nally leading to the molecular-state representation for electronic 
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(fixed-nu clei) wave functions and nonadiabati c nuclear dynam- 
ics (e g. IMott & Massevl 1 1 949t iBatesI Il962t iMacfas & Rienl 
11982k iNikitin & Umanskvlll984l) . Total wave functions of the 
entire collisional system are expressed in terms of products of 
electronic molecular-state wave functions iffi°Hr; R) and nuclear 
wave functions, e.g., radial nuclear wave functions Fi( U )(R)/R 
multiplied by angular nuclear wave functions; R being nuclear 
coordinates. This full quantum treatment requires calculations 
of quantum-chemical data (molecular potentials and nonadia- 
batic couplings) and integration of coupled-channel scattering 
equations. The nuclear motion induces nonadiabatic transitions 
between molecular states. In the asymptotic region, where the 
internuclear distance is large enough, molecular orbitals are de- 
termined by corresponding atomic orbitals, but at short and in- 
termediate distances, where nonadiabatic transitions mainly take 
place, molecular orbitals are usually represented by mixtures of 
quite different atomic orbitals, for example, covalent and ionic 
configurations. 

Perturbation theory can also be applied to low-energy atomic 
collisions under certain conditions. This can be done in the 
adiabatic or diabatic representations. In both representations 
the unperturbed functions are molecular-state electronic wave 
functions, adiabatic or diabatic ones, and perturbation operators 
appear from the nuclear kinetic energy operator. In the adia- 
batic represent ation, the transition probability amplitude reads 
(Belyaev 2007) 



1 f°° 

2 1 « 



F'f(R) 



d_ 
dR 
dFf(R) 
dR 



, dFf 
-Ff(R) ' 



(R) 



dR 



dR, 



(7) 



where (i//'™ l \d / dR\iff'™ l ) r is the nonadiabatic coupling, which 
is responsible for nonadiabatic transitions between adiabatic 
molec ular states / and u due to the nuclear motion; see Belyaev 
(2007) for details. It is clearly seen from Eq.|7]that inelastic tran- 
sitions in low-energy atomic collisions are determined by the nu- 
clear kinetic energy operator (i.e. -ft 2 Vg/2yu), by the mixing of 
electronic molecular-state wave functions expressed by nonadi- 
abatic couplings, and by nuclear wave functions, which substan- 
tially deviate from plane waves. In the dia batic representation, 
the co rresponding amplitude has the form (Niki tin & Umanskvl 
[1981 



i diab 



Ff^WH/uF^WdR , 



(8) 



where H/„ is an off-diagonal matrix element calculated 
with diabatic electronic molecular-state wave functions 
dah i//™y which are obtained from adiabatic quantum-chemical 
data by the diabatisation procedure u nder the condition 
<* a Vr 1 l <5 /dfi| diab C ol >r = dSmithlll969l) . Again, the basis for 
nonadiabatic transitions is the nuclear kinetic energy operator, 
the molecular-state wave functions and nuclear wave functions, 
which should be calculated in this case with diabatic potentials. 

Comparison of Eqs. (0 and ^ for low-energy collisions 
with Eqs. Q and (HJi for high-energy collisions shows that the 
basic physics behind the slow atomic collision theory and the 
first Born approximation, which can be related to the classical 
description, is different. In high-energy collisions the perturba- 
tion theory is based on the unperturbed atomic-state wave func- 
tions and an interaction between a charged impact particle and 
atomic electrons, while in low-energy atomic collisions the per- 
turbation theory is based on the molecular-state wave functions 



and an operator appearing from the nuclear kinetic energy oper- 
ator. If one would use unperturbed electronic atomic-state wave 
functions within perturbation theory for slow atomic collisions, 
a number of new terms related to other interactions would appear 
in addition to the operators mentioned above. 

Thus, it is not surprising that, as we will see, the Drawin for- 
mula and the quantum calculations give substantially different 
results. In some cases, simplifi ed nonadiabatic t ransition mo dels, 
e.g., the Landau -Zener m odel (lLandaul !932a b; Zener 1932J), the 
Demkov ( 19641) . and the INikitin! d 19621) ones, which are based 
on the molecular-state representation, can give reliable inelastic 
cross sections. Nevertheless, the most accurate results are ob- 
tained by means of complete quantum treatments. 

2.3. Quantum collision calculations 

The basic theory of atomic collisi ons has been reasonab l y well 
under stood for some time (e.g. IMott & Massevl 119491: IBatesI 
1 19621 and references therein). However, the detailed calculation 
of accurate quantum-chemical data for quasi-molecules required 
as input to the relevant coupled-channel scattering equations, 
including for excited states, has been a major hurdle. Further, 
while the solution of the scattering equations in the standard adi- 
abatic approach appears straightforward, it encounters difficul- 
ties such as the electron translation prob lem (Bates & McCarroll 
ll958tlBelvaev et al.l[200lt iBelvaevlbo 1 Ot) . 

Over the last 35 years, our understanding of diatomic hy- 
dride molecules and low-energy H atom collisions with atoms, 
particularly alkalis, has increased markedly. Detailed theoretical 
stu dies of alka li hydride molecular structure began with those 
bv ISachs et all dl975l); iMever & Rosmusl (119751) ; lOlson & Liul 
( 19801); IM6 et all dl985l)T and these works showed the impor- 
tance of the ionic configura tion, as had be en expected from 
experimental studies (e.g. Mul likenl 1 19361) . Studies of ion- 
pair production and mutual neutralisation processes in alkali- 



hydrogen collisions fol lowed shortly afte r (Janev & Radulovic 
1 19781; Olson & Kimura|[T985l; lErrea et al.lll986l:lDickinson et alj 
[19991 ICroftet al.ll 19991) 

Though lErrea et all (|l986) briefly looked at excitation of 
sodium by high-energy H atom collisions, the study of excita- 
tion processes came s omewhat la ter. Motivated by Steenbock 
and Holweger's work, iFleck et all (Il99ll) performed an exper- 
imental study of Na(3s) + H — > Na(3/?) + H at low energies (15- 
1500 eV), though due to experimental difficulties, not down to 
near the threshold (2.1 eV), which is the relevant regime for the 
te mperatures of i nteres t. The results were compared with theory 
by iMo & Rieral (1 19921) at energies above 200 eV. Revised ex- 
perimental data, including results down to 10 eV, and detailed 
quantum scatt ering calcul a tions down to the threshold, were 
presented in lBelvaev et al.l (fl999). These studies demonstrated, 
through comparison with the Landau-Zener model and detailed 
quantum scattering calculations, that the experimental data could 
be reasonably explained by nonadiabatic transitions associated 
with avoided crossings, also called pseudocrossings, in the NaH 
molecular potentials. At pseudocrossings, adiabatic potentials 
approach each other and the transition probabilities may become 
large even for slow collisions, despite that the Massey criterion 
based atomic energy level spacin gs predicts the collisions to be 
nearly adiabatic (e.g. lBate "sl[l962, pg. 582 and 611). 

In the particular case of low-lying states of NaH, as discussed 
above, these pseudocrossings arise from interaction of the cova- 
lent Na+H and ionic Na + +H~ configurations, and are thus called 
avoided ionic crossings. The adiabatic potentials showing these 
pseudocrossings are shown in Fig. [2] where we see clearly the 
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influence of the ionic configuration, which strongly perturbs the 
potentials from their behaviour at large internuclear distance. 
Simil ar effects are see n in other quasi -molecules such as LiH 
(e.g. ICroftetaIlll999l) and MgH (e.g. iGuitou et~ai1l201 ll) . As 
shown by diabatic studies, the underlying ionic curve behaves 
roughly as 1 /R and pseudocrossings occur with all neutral states 
below the i onic limit, which show a much weaker dependenc e 
on R (e.g. iBoutalib & Gaded fl992h iBerriche & Gadeal 1 19951) 
Quan tum scatterin g calculations for Li +H ([Belvaev & Barklem 
120031) and Na+H dBelvaev et alll2010h including all states be- 
low the ionic state asymptotic limit have found this to be an 
important mechanism for nonadiabatic tran sitions. Rec e nt cal - 
culations for inelastic Mg+H collisions bv lGuitou et alJ (1201 lh . 
though presently only done for a small number of states, also 
find this result. 

The basic physical interpretation of this mechanism is that, 
in a given A+H quasi-molecule, starting from a covalent state 
with appropriate symmetry where the electrons are associated 
to their respective atoms, at avoided ionic crossings the valence 
electron associated with the atom A has a probability to tunnel 
to the H atom resulting in a predominantly ionic charge distri- 
bution A + +H . Later during the collision at a different avoided 
crossing, there is a probability that the electron may tunnel back 
to a different covalent molecular state leading finally to a dif- 
ferent final state of the atom A, and excitation or deexcitation 
A(nl) + H — > A(n'Z') + H. The electron may also stay with 
the H atom leading to ion-pair production A + H — > A + + H . 
Calculations for Li+H and Na+H mentioned above have gener- 
ally shown that large-valued cross sections for excitation are typ- 
ically small compared with those for ion-pair production from 
certain states. This can be easily understood: a large-valued 
cross section for an excitation process results from passing two 
avoided ionic crossings, one of them (higher-lying) passed with 
a small nonadiabatic transition probability, while for an ion-pair 
production process this region is passed with a large transition 
probability due to small adiabatic energy splittings at highly ex- 
cited crossings. 

We conclude by emphasizing that the mechanism found to 
be important in Li+H, Na+H and Mg+H collisions discussed 
above, that of the avoided ionic crossings, is an example of more 
general phenomena in atomic collisions, wh ich will be discussed 
further below. In He+H inelastic collisions Belvaev et al. (2001) 
couplings arise due to electronic structure properties completely 
unrelated to the ionic configuration. Syste ms such as Li+N a, on 
the other hand, show a similar mechanism (lBelvaevll2010h . 

2.4. Discussion 

Low-energy atomic collisions differ at a fundamental level from 
electron collisions, especially those at high energies. This differ- 
ence arises, naturally, from the presence of a complex projectile 
and thus the effects o f nuclear dyn amics in the atomic collisions. 
As put succinctly bv lBatesi dl962t) : "They are more complicated 
in that a wider variety of reaction paths may be followed: for ex- 
ample, the projectile as well as the target may have structure and 
hence may undergo excitation or ionization; and again, in ad- 
dition to the possibility of electron exchange, the possibility of 
charge transfer arises." Put another way, the effect of the nuclear 
dynamics is that the electronic structure of the atom-atom system 
during the collision can differ significantly from the electronic 
structure of the separated atoms, leading to interactions and cou- 
plings that are very different in character to those found in the 
separated atoms. The phenomenon of the avoided ionic cross- 
ings, where an electron may tunnel from one nucleus to the other 
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Fig. 2. Adiabatic potential energies as a function of internu- 
clear distance R for the lowest ten 'E + states o f the NaH quasi- 
molec ule from pseudo-potential calculations by Dick inson et alJ 
(1999). The atomic states and energies at dissociation are shown 
at the right-hand side of the figure. The thick, dashed (red) line 
shows the 1 /R interaction corresponding to the interaction for 
the pure ionic configuration Na + +H" at long range. 



to form an ionic electronic structure leading to strong coupling 
between adiabatic states, is a particular example of the effects of 
nuclear dynamics in atomic collisions. In contrast, electron colli- 
sions with atoms have no such possibility to change the structure 
of the target atom in the same way. 

Thus, we see there is little in common between the physics 
of the Drawin formula and the physics of low-energy atomic 
collisions, especially H atom collisions, as we now understand 
them. The Drawin formula is based on a classical model and we 
see that the fundamental mechanism involved in H atom colli- 
sions for the cases studied relates to electron tunnelling, which 
is quantum mechanical in nature. The Drawin formula extends 
a theory for electrons, neglecting the atomic nuclei other than 
in scaling the collision energy. Thus, from a physics point of 
view it seems the grounds for expectations that Drawin formula 
can provide good estimates are rather weak, despite the fact that 
Thomson formula on which it is based is capable of providing 
order-of-magnitude estimates in some situations for ionization 
by electron impact. 
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3. Comparison of results 

Detailed quantum scattering calculations of cross sections, and 
corresponding rate coefficients, have been done for excitation 
processes between all states below the ionic lim it for Li+H 
(Belv aev &Barklern] 12003b IBarklem et al l 120031) and Na+H 
(Bel vaev et alj l201(t IBarklem et all |2Q10|> . Recently, calcula- 
tions fo r the three lowest s tates of Mg+H have also been per- 
formed (Guitou et aD l201 ll) . In this section we compare the re- 
sults with the predictions of the Drawin formula. This compari- 
son will be done on two levels. Firstly, comparison of the cross 
sections, which best elucidates the physics. Secondly, compari- 
son of the rate coefficients, which is of most interest for astro- 
physical applications, and also allows us to focus on more gen- 
eral integrated properties. 
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3.1. Cross sections 

In this section we compare the Drawin cross sections 
with other theoretical and e xperimental results. Note, since 
ISteenbock & Holwegerl d!984l) give only a formula for the rate 
coefficient, and due to above mentioned problems with Drawin's 
deriv ation, the Drawi n cross section formula used here is taken 
from [Lambertl (Il993l Eqs. A2 and A8). As shown there, these 
expressions give a formula for the rate coefficient differing 
from that of ISteenbock & Holwegerl (|l984) by only a factor of 
twa/Owa + '«hX which is of order unity. 

Fig. [3] compares cross sections at low energy for Na(3s) + 
H — » Na(3p) + H, the single case where experimental data is 
available. We see that quantum scattering calculations and the 
Landau-Zener model results agree quite well with experiment. 
Near the threshold, where there is no experimental data, the 
quantum scattering calculations show substantial differences de- 
pending on which quantum-chemical data are used. This high- 
lights the sensitivity of the near-threshold cross sections to the 
uncertainties in the quantum-chemistry data, and gives an es- 
timate of the uncertainties in the calculated cross sections. In 
this case the uncertainties seem to be around one or two orders 
of magnitude. However, it should be noted that the uncertain- 
ties vary strongly from transition to transition, and the transi- 
tions with the largest cross sections have the smallest uncer- 
tainties, as low as a factor of 2. Fortunately, these are the most 
important from an a strophysical perspective (see discussion in 
IBarklem et af]|2010h . 

In Figs. |4] and [5] theoretical excitation cross sections for in- 
elastic Li+H and Na+H collisions at low energy are compared 
for a number of transitions. A major problem of the Drawin 
formula is immediately evident here: it predicts zero cross sec- 
tions for optically forbidden transitions (or cross sections much 
smaller than those for optically allowed transitions if one uses 
the very small /-values in some cases). The Drawin formula de- 
pends on the absorption oscillator strength, while collisional pro- 
cesses do not follow the same selection rules, and no such depen- 
dence is seen in the quantum mechanical data. For the optically 
allowed transitions the Drawin cross section is typically greater 
than the quantum scattering results by several orders of mag- 
nitude. The Drawin formula gives a reasonable description of 
the cross section behaviour with collision energy in many cases, 
though in a few, particularly among transitions involving excited 
initial states, the quantum cross sections have a much flatter be- 
haviour with collision energy. 



Fig. 3. Comparison of cross sections for the process Na(3s) + 
H — > Na(3p) + H. The dotted line shows the Drawin cross sec- 
tion and the dashed line the Landau-Zener cross section. The 
full lines are two quantum scattering calculations using differ- 
ent input quantum-chemical data (MR DCI, which is larges t near 
the threshold, and pseudopotential; see Bel vaev et al.l(l2010t) for 
details ) . The circles show th e experimental data of Fleck et alj 
(ll99lh : lBelvaev et all (Q999) with lcr error bars. 



3.2. Rate coefficients 

In Fig. |6l rate coefficients for excitation at 6000 K by in- 
elastic Li+H, Na+H and Mg+H collisions, found by integra- 
tion of the cross sections over a Maxwellian velocity distri- 
bution, are shown and compared. In addition to the quantum 
mechanical calculations and the Drawin formula results, we 
also calculat e rate coefficient s from the free electron model 
described by iKaulakvsl d!99 ll) . Cross sections were calculated 
according to Eq. 18 in Kaulakys, using non-hydrogenic wave 
functions in momentum space calcula ted using the methods of 
iHoang Binh & Van Regemorterl(ll997l) . This model is only valid 
for transitions in Rydberg atoms, and thus is not valid for the 
low-lying states being examined here; however, we make a com- 
parison in any case to see if there is any possibility it may give 
reasonable estimates for transitions among low-lying states. The 
rate coefficients are plotted against the transition energy, AE, 
since this a main parameter in determining the rate coefficients, 
as seen in the figure, due to the strong dependence of the in- 
tegral on the threshold (e.g., see Eq. Q] where the exponential 
dependence on AE = E^ is easily seen). This allows us to sep- 
arate this effect from other more intrinsic secondary effects re- 
lated to the physics of the collisions. The quantum mechanical 
results correlate strongly with AE, though with a significant scat- 
ter of several orders of magnitude around the mean relation. One 
sees, particularly in the case of Na+H, clear secondary trends 
among transitions with the same initial state, e.g. 3p in Na+H. 
These trends result from regular behaviour in the properties of 
the quasi-molecule at avoided crossings (potential splitting, ra- 
dial coupling), leading to trends in the transition probabilities. 

Comparison with the Drawin formula immediately shows 
the problem already discussed: lack of predictions for forbidden 
transitions. The Drawin formula data also show a strong corre- 
lation with AE, as would be expected since for the case of a 
given element and temperature, the formula has only two pa- 
rameters, AE = E^, and /;„. Thus, all deviation from the main 
trend with AE is a result of the /-value and those transitions de- 
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Fig. 6. Rate coefficients at 6000 K from various calculations for excitation transitions between the lowest 9 states of Li and Na due to 
H atom collisions, plotted against the transition energy. The panels in the left column are for Li+H, and in the right column Na+H. 
In each panel, different symbols are used to denote the initial state of the trans ition following the key g iven in the upper panels. The 
upper panels show the results from detailed quantum scattering calculations (Bark lem et al.l l2003, 2010), noting there are no results 
for Li n — 4 initial states as the data for these transitions were deemed unreliable. The dotted line shows a linear fit to the these data, 
which is repeated in the lower two panels to aid comparison. The upper panel in the right column also shows the limited data now 
available for Mg+H calculated from cross sections given in iGuitou et al.l (1201 lb . Note that the data for the two transitions to the 
3p 'P° state should be interpreted as upper limits; however, they are probably not more than an order of magnitude or two lower. No 
data for Mg+H are plotted in the other panels. The middle panels show the results of the Drawin formula, and the insets show the 
ratio with the quantum scattering calculations. Forbidden transitions, where the Drawin form ula is not applica ble, are shown in the 
bottom of the panels. The lower panels show the results of the Kaulakys free electron model (Kaulakv3 ll99lh . and the insets show 
the ratio with the quantum scattering calculations. 
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Fig. 4. Co mparison of quantum scatt ering cross sections for 
Li+H from lBelvaev & BarklerrJ (I2003I) (full lines) with those of 
the Drawin formula (dotted lines) for the 10 transitions between 
the 4 lowest states of Li. 



Fig. 5. Compa rison of quantum scattering cross sections for 
Na+H from iBelvaev et alj (|2010) (full lines) with those of the 
Drawin formula (dotted lines) for the 10 transitions between the 
4 lowest states of Na. 



viating significantly from the trend are those with very high or 
low /-values. These secondary trends are seen not to correlate 
with those seen in the quantum mechanical data, which is ex- 
pected as the physical mechanisms are not related. In general, it 
is seen that the rate coefficients differ from the quantum ones by 
up to eight orders of magnitude with a general tendency to be 
significantly larger. There is also a huge scatter spanning around 
eight orders of magnitude. 

Though the Kaulakys free electron model has the advan- 
tage of providing data for all transitions, it compares no better 
with the detailed quantum scattering results in general, except 
perhaps for transitions with small transition energies AE where 
agreement is within a few orders of magnitude. Comparison 
shows a clear tendency for increasingly poor agreement with in- 
creasing transition energy, with significant scatter. Thus, there 
is no suggestion that the Kaulakys free electron model could 
give reasonable estimates for transitions among low-lying states. 



It should, however, provide reasonable estimates for transitions 
among Rydberg states above the ionic limit. 

4. Conclusions 

The Drawin formula does not contain the essential physics be- 
hind direct excitation by H atom collisions at low energies, 
which is quantum mechanical in character. It compares poorly 
with the results of the available full quantum scattering calcu- 
lations based on detailed quantum-chemical modelling of the 
quasi-molecule. The estimated uncertainties in the rate coeffi- 
cients from the quantum mechan ical data are not su fficient to 
explain the differences (see, e.g., Bar klem et al.ll2010l) . In gen- 
eral, the Drawin formula significantly overestimates the rates by 
amounts that vary markedly from transition to transition. Thus, 
we conclude that the Drawin formula does not provide order-of- 
magnitude estimates for excitation by H atom collisions. We can 
even more generally conclude that formulae based on the Drawin 
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expressions are also unlikely to be useful for excitation by other 
atoms. 

The Drawin formula is often used in non-LTE stellar atmo- 
sphere applications in a semi-empirical fashion with a scaling 
factor, often denoted Sh- The problem of forbidden transitions 
is sometimes avoided with some arbitrary solution to allow data 
to be estimated, for example a n arbitrary constant /-value or 
scaled /-value (e.g. ISteenbocklll985l) . The value of Sh is var- 
ied either to assess the possible influence of H atom collisions 
or to calibrate it against some observational criteria. This proce- 
dure, in the absence of a better alternative, is not unreasonable 
in view of the fact that AE has a strong influence on the rate 
coefficient, and this influence resulting from the velocity distri- 
bution is accounted for by the Drawin formula. However, as we 
have seen, there is a large scatter around this basic behaviour of 
the rate coefficients, which stems from the physics of the atomic 
collision, and is not reproduced by the Drawin formula. This 
means such a procedure will not give good estimates of relative 
rates for collisional transitions, and could lead to modelling er- 
rors if more than one transition is of importance. Furthermore, 
if other atomic collision processes are present but not modelled, 
effects caused by other processes may be incorrectly attributed 
to H atom collisions. For example, in modelling of Li and Na 
non-LTE line formation it has been shown that the ion-pair pro- 
duction processes mentioned in Sect. 12.31 often called charge 
exchange or charge transfer, are far more import ant than direct 
excitation dBarklem et alj|2003t iLiiid" et alj|20l"lb . Deficiencies 
in other aspects of the spectrum modelling might also be masked 
by such an approach. 

It is clear that there is a need for quantum mechanical es- 
timates for direct excitation and ion-pair production processes. 
The most accurate calculations are naturally provided by the full 
quantum treatment, i.e. detailed quantum-chemical data calcu- 
lation with detailed quantum scattering calculations. However, 
in order to provide estimates for a wide range of elements of 
astrophysical interest, simplified models, based on the physical 
understanding provided by the detailed calculations, would be 
required. The reasonable success of the Landau-Zener model, as 
seen in Fig. [3] suggests this model, when coupled with a method 
for producing quantum-chemical data, could lie at the basis of 
attempts to produce such simplified models. However, we em- 
phasise that more experience from the full quantum calculations 
is complementary and still needed. 
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